Identification of the bulk pairing symmetry in high-temperature superconductors: 
Evidence for an extended s-wave with eight line nodes 
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we identify the intrinsic bulk pairing symmetry for both electron and hole-doped cuprates from the 
existing bulk- and nearly bulk-sensitive experimental results such as magnetic penetration depth, 
Raman scattering, single-particle tunneling, Andreev reflection, nonlinear Meissner effect, neu- 
tron scattering, thermal conductivity, specific heat, and angle-resolved photoemission spectroscopy. 
These experiments consistently show that the dominant bulk pairing symmetry in hole-doped 
cuprates is of extended s-wave with eight line nodes, and of anisotropic s-wave in electron-doped 
cuprates. The proposed pairing symmetries do not contradict some surface- and phase-sensitive 
experiments which show a predominant d-wave pairing symmetry at the degraded surfaces. We also 
quantitatively explain the phase-sensitive experiments along the c-axis for both Bi2Sr2CaCu2 08+i/ 
and YBa2Cu307_ !/ . 



I. INTRODUCTION 

An unambiguous determination of the symmetry of 
the order parameter (pair wavefunction) in cuprates is 
crucial to the understanding of the pairing mechanism 
of high-temperature superconductivity. In recent years, 
many experiments have been designed to test the or- 
der parameter (OP) symmetry in the cuprate supercon- 
ductors. However, contradictory conclusions have been 
drawn from different experimental techniques [|l] [l6| , 
which can be classified into being bulk-sensitive and 
surface-sensitive. For example, the magnetic penetra- 
tion depth measurements and polarized Raman scat- 
tering experiments are bulk-sensitive. Angle-resolved 
photoemission spectroscopy (ARPES) is essentially a 
surface-sensitive technique. However, the ARPES data 
for Bi 2 Sr 2 CaCu 2 0s+ a (BSCCO) should nearly reflect the 
bulk properties since the cleaved top surface contains 
an inactive Bi-0 layer, and the superconducting coher- 
ent length along the c-axis is very short. The single- 
particle tunneling experiments can probe the bulk elec- 
tronic density of states when the mean free path is far 
larger than the thickness of the degraded surface layer 
p7[ . Therefore, the single-particle tunneling experiments 
along the Cu0 2 planes are almost bulk-sensitive due to 
a large in-plane mean free path (> 100 A). In contrast, 
the phase-sensitive experiments based on the Josephson 
tunneling are rather surface sensitive (since pair tunnel- 
ing is limited by the coherence length which is rather 
short in cuprates), so that they might not probe the 
intrinsic bulk superconducting state if the surfaces are 
strongly degraded. In this case, the observed product of 
the critical current times the junction normal-state re- 
sistance (I c Rn) will be very small compared with the 
Ambegaokar-Baratoff limit. Then the OP symmetry at 
surfaces may be different from the one in the bulk p8|. 
Therefore, the surface- and phase-sensitive experiments 
do not necessarily provide an acid test for the intrinsic 
bulk OP symmetry. 



Here, we identify the intrinsic bulk pairing symmetry 
for both electron and hole-doped cuprates from the exist- 
ing bulk- and nearly bulk-sensitive experimental results 
such as magnetic penetration depth, Raman scattering, 
single-particle tunneling, Andreev reflection, nonlinear 
Meissner effect, neutron scattering, thermal conductivity, 
specific heat, and ARPES. These experiments consis- 
tently show that the dominant bulk pairing symme- 
try in hole-doped cuprates is of extended s-wave with 
eight line nodes, and of anisotropic s-wave in electron- 
doped cuprates. The proposed pairing symmetries do 
not contradict some surface- and phase-sensitive exper- 
iments which show a d-wave pairing symmetry at the 
degraded surfaces. The extended s-wave pairing sym- 
metry deduced from the bulk-sensitive experiments is 
also in quantitative agreeement with the well-designed 
phase-sensitive experiments along the c-axis for both 
Bi 2 Sr 2 CaCu 2 0s+ a and YBa 2 Cu307_ y . 



II. THE PAIRING SYMMETRY IN HOLE-DOPED 
CUPRATES 

A. The pairing symmetry in Bi2Sr2CaCu20s+y 

We first examine the high-resolution ARPES data ob- 
tained for Bi 2 Sr 2 CaCu 2 08+ y crystals 14 From the 



ARPES data, one can determine the angle dependence of 
the superconducting gap with a resolution as high as ± 2 
meV |l(| . Due to the complication arising from a possible 
superlattice contribution in the X quadrant, we only use 
the data obtained for the Y quadrant to extract the gap 
function. In Fig. ^, we show the angle dependence of the 
superconducting gap A(0) in the Y quadrant for slightly 
overdoped and heavily overdoped BSCCO single crystals. 
The data were taken from Ref . |l4],[l6| . Here 9 is the an- 
gle measured from the Cu-0 bonding direction. For the 
slightly overdoped sample (Fig. |l|a), the gap at 9 
= 45° (diagonal direction) is very small (3.5±2.5 meV), 
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and the gap symmetry could be consistent with a d-wave 
symmetry, i.e., A(6) = A cos 26. On the other hand, 
the gap along the diagonal direction (r~Y) for the heav- 
ily overdoped sample (Fig. [l]b) is not small (9±2 meV), 
which is obviously not consistent with the d-wave pairing 
symmetry. A similar evolution of the gap function with 
the doping has been observed by the bulk-sensitive po- 
larized Raman scattering 0, which also shows that the 
difference in the magnitudes of the gaps along the Cu-0 
bonding direction and the diagonals becomes smaller and 
smaller towards overdoping. 
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FIG. 1. The angle dependence of the superconducting gap 
A(0) in the Y quadrant for BiaSraCaCuaOs+y (BSCCO) crys- 
tals: (a) slightly overdoped sample with T c = 87 K; (b) heavily 
overdoped sample with T c — 60 K. The magnitudes of the gap 
were extracted from ARPES data Here 6 is the angle 

measured from the Cu-0 bonding direction. 

The question is: what functional form of A(6) can 
fit the angle dependence of the gap shown in Fig. [I]? 
In general, the gap can be expressed as A(0) = A s + 



A d cos 26 + A g cos 46 
has 



In the case of A^ ~ 0, one 



A(0) = A(cos40 + s), 



(1) 



where s is the parameter reflecting the isotropic s-wave 
component. This gap function has eight line nodes for 
s < 1, while there are no nodes for s > 1. The gap 
function (Eq. |l|) is also called extended s-wave (denoted 
by s*-wave). The polarized Raman data for an optimally 
doped HgBa2CaCu206+j/ are in good agreement with the 
s*-wave gap function M. If we take the absolute value of 



A (6>), then 

|A(0)| = |A(cos40 + s)|. (2) 

We fit the data of Fig. |l| by Eq. ||. It is remarkable that 
the fits are rather good. This indicates that the ARPES 
data may be consistent with the extended s-wave sym- 
metry. The ARPES specified maximum gap Am at 6 = 
for the slightly overdoped sample is 36±3 meV, which 
is much larger than the value (~28 meV) determined 
from break junction spectra ^J. On the other hand, 
the ARPES determined A M value (15±2 meV) for the 
heavily overdoped sample with T c = 60 K is very close 
to the value (18±2 meV) inferred from a break junction 
spectrum of a similar crystal with T c = 62 K p2[ . The 
discrepancy in the former case may be due to the fact 
that the doping level in the top layer where the ARPES 
probes could be slightly lower than in the bulk (i.e., the 
top Cu02 layer might be slightly underdoped). Thus, the 
ARPES experiments on the BSCCO single crystals are 
nearly bulk sensitive, in contrast to the ARPES exper- 
iments on other cuprates, which are essentially surface 
sensitive. 

If the proposed gap functions (Eq. [I] and Eq. ||) are 
indeed relevant, they should be also consistent with other 
bulk-sensitive experimental results such as the in-plane 
magnetic penetration depth X ab (T). Since there are eight 
line nodes in the proposed gap function, the change of the 
in-plane penetration depth at low temperatures should be 
proportional to T. Following the procedure in Ref. [p3| , 
we can readily show that the slope 



d\ ab (T)/dT = [A o6 (0) ln2/A M Ml + s)/(l-s). (3) 

Compared with the (i-wave symmetry, the magnitude 
of the slope dX ab (T)/dT is enhanced by a factor of 
y(T~+ s)/(l — s). In terms of Am and Ad, we find that 
s = (A M - A D )/(A M + Ad) and A = (A M + A D )/2. 
Then, Eq. || can be rewritten as 

d\ ab {T)/dT = X ab (0) In2/ V / A M A^. (4) 

It is interesting to see that dX ab (T) /[X a b{0)dT] is in- 
versely proportional to yAyAo, namely, the geometric 
average of Am and A^. 

The single-particle tunneling spectroscopy can probe 
the superconducting density of states (DOS) with fine 
energy resolution and considerable directionality. For 
an isotropic s-wave superconductor, the characteristic 
dl/dV vs V curve in the point-contact SIN tunneling 
junctions exhibits a step- like peak at a voltage V p = A/e. 
For an anisotropic gap function A(6), the directional de- 
pendence of the tunneling differential conduction is given 

bypl 



2^ 



^« / m6-M[ 



eV ~iT 



\N{8)d8. 
(5) 
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Here N(9) represents the anisotropy of the band disper- 
sion; r is the life-time broadening parameter of an elec- 
tron; p(0 — 60) is the angle dependence of the tunneling 
probability which decays exponentially as p(0 — 9o) = 
exp[-/?sin 2 (0-0 o )] (60 is the angle of the tunneling bar- 
rier direction); the parameter (3 decreases with decreas- 
ing barrier resistance Rn- For simplicity, we assume a 
cylindrical Fermi surface, so that both N(0) and (3 are 
independent of the angle. This will not change the ba- 
sic features of the dl/dV curve. In Fig. || we show the 
numerically calculated results of the renormalized dl/dV 
for a gap function of A(9) = A(cos4# + s) with A = 24 
meV and s = 0.25. One can readily show that the max- 
imum gap is Am = (1 + s)A = 30 meV at 9 = 0, and 
the gap along the diagonal directions is Ad = (1 — s)A 
= 18 meV. From Fig. ||, one can see that either two or 
four peak features appear clearly in the dl/dV curves, 
depending on the tunneling barrier direction and/or the 
j3 value. For a small (3 value (corresponding to a small 
barrier resistance), four peak features are well defined 
(see curve A). The peak positions are located at eV = 
±Am and ±Ae>. Therefore, from the peak positions, we 
can determine Am and A^. 
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FIG. 2. Numerically calculated curves of the renormalized 
dl/dV for a gap function of A(t9) = A(cos4# + s) with A = 
24 meV and s = 0.25. The four curves correspond to different 
values of the parameters F, j3, and o , which are indicated in 
the figure. The curves A, B and C are vertically shifted up 
by 3, 2 and 1, respectively. 

In Fig. |^, we plot the normalized dl/dV curve at 14 
K for an SIS break junction on a BSCCO crystal which 
is slightly overdoped (T c = 90 K) |2jJ. The junction has 
a very low barrier resistance (~200 Q) j25| , indicating 
a small (3 value. It is remarkable that there are four 



well-defined peak features in the spectrum, which resem- 
ble curve A in Fig. ^. The pronounced zero bais peak 
arises from Josephson tunneling [^T],^5j . When the bar- 
rier resistance is above 2 kfi, the inner gap features dis- 
appear p5|] , in agreement with curve B in Fig. § We 
would like to mention that, for SIS break junctions, the 
peak positions are located at eV = ±2 Am and ±2A£>. 
From the spectra, we obtain Am = 26±0.5 meV, and 
An = 9.5±0.5 meV. The Am value obtained from the 
break junction spectrum is the same as that found from 
the c-axis intrinsic tunneling junctions made of the in- 
sulating Bi-0 layers |26j. From the Am and Ab values, 
we deduce a gap function A(0) = A(cos46* + s) with 
A = 17.75 meV and s — 0.46. With this gap function 
and A a6 (0) = 2690±150 A§,we calculate from Eq. 4 
that d\ ab (T)/dT = 10.2±0.6 A/K, in excellent agree- 
ment with the measured values (10.2±0.2 A/K) gg. 
Similarly, the earlier break junction spectra for an over- 
doped BSCCO with T c — 86 K also indicate double gap 
features at Am = 24±2 meV and at Ad = 12±1 meV 
(Ref. H^]). The tunneling spectra are in good agreement 
with ARPES data for an overdoped BSCCO with T c = 
83 K @. The ARPES experiment clearly showed that 
Am = 20±2 meV and A D = 12±2 meV @. Moreover, 
the inner gap features also appear in SIS break junc- 
tion spectra of a heavily overdoped crystal with T c — 
62 K, corresponding to A D = 7.5-9.0 meV (Ref. §|||). 
The magnitude of Ac is in excellent agreement with that 
found from the ARPES experiment (see Fig. |l|b) . 
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FIG. 3. Normalized dl/dV curves at 14 K for the SIS break 
junctions on a slightly overdoped BSCCO crystal. The spec- 
tra were taken from Ref. pB| . 

We would like to point out that the values of Am de- 
termined from Raman spectrum of B\ g symmetry may be 
overestimated due to the fact that the extended van Hove 
singularity is slightly below the Fermi level. In this case, 
the spectra would show double peaks at Raman shifts of 
2 Am and 2-\/A| f + £% H , where £ vH is the energy posi- 
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tion of the van Hove singularity below the Fermi level. 
When £ v h « Am, one can only see a single broad peak 
slightly below 2y/A 2 M +g H . 



B. The pairing symmetry in YBa2Cu307_ y 

Evidence for an extended s-wave pairing symmetry in 
YBa 2 Cu 3 07_ a (YBCO) also comes from single-particle 
tunneling spectra. Fig. |^ shows scanning tunneling spec- 
trum for a slightly overdoped YBa2Cu30y_ y crystal [ p0[ . 
Four peak features appear in this spectrum that is sim- 
ilar to curve D in Fig. |^. From the peak positions, we 
obtain A M = 30±2 meV, and A D = 19±1 meV. The 
size of Am — 30 meV is consistent with a break junction 
spectrum JT^ ] , and a scanning tunneling spectrum along 
the a-axis direction H. A gap feature with Ac = 19 
meV was also seen in a scanning tunneling spectrum Q 
that is very similar to curve C in Fig. ^[ 
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FIG. 4. Scanning tunneling spectrum for a slightly over- 
doped YBa2Cu3C>7 (YBCO) crystal. The spectrum was taken 
fromRef. @. 

Now we discuss the Andreev reflection. Since there is 
sign change about its nodal directions in our extended s- 
wave order parameter, the Andreev-bound surface states 
can be formed. This will lead to a zero-bias conduc- 
tion peak (ZBCP) if tunneling is nearly along one of the 
nodal directions and the bare Fermi velocities between 
the cuprates and normal metals (e.g., Ag and Au) are 
well matched. For hole doped cuprates, the bare Fermi 
velocity v b F strongly depend on the angle 9, that is, v F 
is small along the bonding direction, and large along the 
diagonal directions. This implies that the observation of 
the Andreev reflection is difficult for tunneling along the 
bonding direction since the value of v F along this direc- 
tion is small compared with that of Au or Ag. Due to 
the strong anisotropy of v F in cuprates, the Andreev re- 
flection mainly probes the gap feature at eV — Ad- If 
tunneling is along one of the diagonal directions, and the 
angle between the nodel and diagonal directions is far 
larger than the half tunneling angle (depending on /?), 
one can see an s-wave like gap approximately equal to 



Ad in the Andreev reflection spectra. Indeed an s-wave 
like gap feature at eV ~ 20 meV has been observed in 
the Andreev reflection spectra of several YBCO crystals 
with T c = 90 K |HJ. We would like to mention that, 
in general, the double gap features should also appear in 
the Andreev reflection spectra when the j3 value is small 
and v b F does not have a significant anisotropy. 

The tunneling data of YBCO (Fig. @) are thus consis- 
tent with a gap function A(0) — A(cos4t9 + s) with A = 
24.5 meV and s = 0.225. This gap function is in quan- 
titative agreement with the a-axis A a (T) data (which re- 
flect magnetic screening in Cu02 planes) for a fully oxy- 
genated YBCO crystal p2[ . From Eq. [|, we calculate 
d\ a {T)/dT = 4.0 A/K using A Q (0) = 1600 A (Ref. |§), 
Ad = 19 meV, and Am = 30 meV. We will get the same 
value of dX a {T) j dT if we use Ad — 21 meV, and Am 
= 27 meV. For a d-wave gap function A(8) = Am cos 29 
with A M = 30 meV, the calculated d\ a {T)/dT = 3.2 
A/K. The measured value of d\ a {T)/dT is 4 A/K ||. 
It is evident that the extended s-wave gap function is in 
much better agreement with experiment than the d-wave 
gap function. 
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FIG. 5. Temperature dependence of the a-axis A„(0)/A„(T) 
for a very high-quality YBa2Cua07 (YBCO) crystal with T c 
= 88.7 K. The data were taken from Ref. |32|. The solid line 
is the calculated curve for the s*-wave gap function deduced 
from the tunneling spectrum in Fig. The dash line is the 
calculated curve for a d-wave gap function with Am = 30 
meV. 

Now we calculate the temperature dependence of 
A^ b (0)/A^ b (T) for the s*-wave gap function. For a cylin- 
drical Fermi surface || 



&(0) 



1 + (1/tt) 



2tt 



Of 

d9dejL. (6) 



Here E = ^e 2 + A 2 ((9,T); / is the Fermi-Dirac distri- 
bution funct ion; A(6>, T) = A(T)(cos46» + s); A(T) = 
Atanh(2.2^/T/r c - 1) (Ref. |g§). In Fig. |, we com- 
pare the experimental data for YBCO (open circles) |3l] 
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and the numerically calculated result (solid line) for the 
above deduced gap function A(0) = 24.5(cos46> + 0.225) 
meV. It is remarkable that the data are in quantitative 
agreement with the calculated result without any fitting 
parameters. The dash line is the calculated result for a 
d-wave gap function A(0) = Am cos 29 with Am = 30 
meV. It is clear that the agreement between the data 
and the calculated curve is poor for the d-wave symme- 
try. It is worthy to note that the temperature depen- 
dence of A^ b (0)/A^ b (T) is mainly determined by the gap 
function, so the shape of the Fermi surface has little effect 
on A 2 fc (0)/A 2 fc (T). 
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FIG. 6. Sine Fourier amplitudes of the transverse mag- 
netization rriT in the Meissner state for a high-quality 
YBa2Cu3C>7 (YBCO) crystal. The data were taken from 
Ref. 0|. The solid line is the predicted sine Fourier amplitude 
at 27r/4 for the s*-wave gap function deduced from the tun- 
neling spectrum in Fig. [ijand the a-axis A a (T) data in Fig. |[ 
The calculated amplitudes at 27r/2, 2n/3, 27r/4, and 2n/5 are 
similar while the ones at other periods are much smaller. The 
dash line is the predicted sine Fourier amplitude at 27r/4 for 
a d-wave gap function. 

The gap function of YBCO deduced from the tunneling 
and the A a (T) data is also consistent with the measured 
transverse magnetization mp in the Meissner state Q , as 
plotted in Fig. || This bulk sensitive experiment shows 
a very small sine fourfold component of the transverse 
magnetization, which is at least 4 times smaller than the 
predicted value from the d-wave symmetry. This indi- 
cates that the dominant pairing symmetry is not the d- 
wave. Using the formulas reported in Ref. pjfl , we can 
calculate the sine components of the transverse magneti- 
zation for the s*-wave gap function deduced above. We 
find that the sine fourfold component for the s*-wave OP 
is a factor of 8.9 smaller than for the pure d-wave OP. 
The predicted sine Fourier amplitude at period 27r/4 is 
indicated by a horizontal solid line in Fig. |g. The cal- 
culated amplitudes at 27r/2, 2ir/3, 2ir/4, and 2ir/5 are 
similar while the ones at other periods are much smaller. 



It is clear that the predicted amplitudes at all the pe- 
riods are below the noise level which is about 5xl0~ 10 
emu Therefore, the very small nonlinear Meissner 
effect observed in the overdoped YBCO is in agreement 
with the s*-wave OP or with a nodeless OP (J] rather 
than with the d-wave OP. A nodeless OP symmetry is in 
contradiction with the observed linear T dependence of 
the thermal conductivity down to a very low temperature 
(50 mK) §|. 

In addition, we further show that the s*-wave gap func- 
tion is in quantitative agreement with the low tempera- 
ture thermal conductivity, specific heat, and surface An- 
dreev bound states. By replacing Am by vAjifAfl in 
the equations for the low-temperature electronic thermal 
conductivity K e i and specific heat C e i for the d-wave gap 
function in the clean limit p5| , we obtain the following 
equations for an s*-wave gap function: 



Kel 

T 



igVpkj- 



6dy/A M A D 



and 



Cei_ = 9C(3)fc|fc f 
T 2 -Khvpdy/A M A D 



(7) 



(8) 



Here vf and kp are the Fermi velocity and momentum 
along the nodal directions, respectively; d is the average 
interlayer distance; £(3) = 1.20. One should note that im- 
purity scattering tends to suppress the values of both 5fr 
and The Fermi velocity along the nodal directions 
has recently been obtained for YBCO from the studies of 
surface Andreev bound states The deduced Fermi 

velocity vf is (I.2±0.2)x 10 5 m/s, which is a factor of 
2 smaller than the measured Fermi velocity along the 
diagonal directions from the ARPES data of BSCCO 
p7| . This suggests that the nodal directions might be 
far away from the diagonal directions. For the s*-wave 
gap function deduced above for overdoped YBCO, the 
nodal directions are about 19° away from the diagonal 
directions (i.e., at 9 = 26°). Indeed, from the ARPES 
data of BSCCO |57j , one can clearly see that the Fermi 
velocity at 9 — 26° is smaller than that at 9 = 45° by a 
factor of about 2. Substituting vf = 1.2xl0 5 m/s, kp 
= 0.7 A^ 1 d = 5.85 A, A M = 30 meV and A D 

= 19 meV into Eq. and Eq. g, we obtain ^ = 0.12 
mW/K 2 cm and = 0.24 mJ/moleK 3 . The calculated 
values are in excellent agreement with the measured val- 
ues: = O.I4±0.03 mW/K 2 cm (Ref. ||) and 
0.20±0.05 mJ/moleK 3 (Ref. 

Moreover, thermal conductivity of YBCO as a func- 
tion of angle of an inplanc magnetic field relative to the 
crystal axes has been studied both theoretically and ex- 
perimentally |3^,|f|. A theoretical calculation for the 
angular dependence of the magnetothermal conductiv- 
ity [ |39"[ shows that an extended s-wave gap produces a 
more symmetric angular variation than a d-wave gap. It 
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appears that both sets of experimental data |3|,|lJ are 
more consistent with an extended s-wave gap than a d- 
wave gap. 



C. The pairing symmetry in La2-a;Sra;Cu04 

The polarized Raman scattering data |y] for nearly 
optimally-doped L^-zSr^CuC^ (LSCO) with T c = 37 
K yield 2Am /ksTc = 7.7. From the measured value of 
d\ ab {T)/[X ab (0)dT} for the o ptimally -doped LSCO @, 
one can readily calculate 2^/ AuAo/kBTc = 4.2 using 
Eq. |. Then we get 2A D /k B T c = 2.3, i.e., A D = 3.8 
meV. This value is in good agreement with the Andreev 
reflection spectrum of optimally-doped LSCO Q, which 
shows the s-wave like gap feature at eV ~ 3.5 meV. 
Therefore, three independent bulk-sensitive experiments 
on the optimally- doped LSCO consistently suggest a gap 
function: A(0) = 8.1(cos40 + 0.53) meV with A D = 3.8 
meV and Am = 12.5 meV. 

Now we can quantitatively explain the neutron scat- 
tering experiment on an optimally-doped LSCO single 
crystal jl4| . The experiment shows that low energy mag- 
netic excitations are peaked at the quartet of wavevectors 
(0.5±0.135, 0.5) and (0.5, 0.5±0.135) in the normal state, 
and a spin gap with energy of about 6.7 meV appears in 
the low-temperature superconducting state. The magni- 
tude of the spin gap should be equal to twice the super- 
conducting gap along the incommensurate wave vectors 
(i.e., at = 39°) From the gap function deduced 

above, we calculate 2A(39°) = 6.2 meV, in remarkably 
good agreement with experiment. Moreover, it was also 
found Q that the spin gap at = 45° is 6±2 meV, which 
is consistent with 2Ajj = 7.6 meV within the experimen- 
tal uncertainty. Obviously, the <i-wave gap function is 
incompatible with the large spin gap observed along the 
diagonal direction. The neutron data might be also con- 
sistent with an isotropic spin gap, as suggested by Lake 
et al. fjij] . However, the isotropic spin gap is incompati- 
ble with the T 3 dependence of the spin-lattice relaxation 
rate observed in hole-doped cuprates. Only with the s*- 
wave gap function for LSCO, one can quantitatively ex- 
plain the neutron experiment, Raman scattering, mag- 
netic penetration depth, Andreev reflection, and mag- 
netic resonances. 



III. THE PAIRING SYMMETRY IN 
ELECTRON-DOPED CUPRATES 

The recent measurements of A a &(T) in an electron- 
doped Pri. 85 Ceo.i5Cu04_ y (PCCO) reveal contradictory 
results 46,4?]]. In a high-quality PCCO thin film with the 
lowest residual resistivity and the highest T c , the temper- 
ature dependence of [X ab (T) — X ab (0)]/ X ab (0) is consistent 
with an s-wave pairing symmetry with a reduced energy 



gap 2A(0)/fc B T c = 2.9 J4§. On the other hand, the low- 
temperature X ab (T) in less ideal PCCO single crystals ex- 
hibits a power-law temperature dependence, as expected 
from a dirty d-wave superconductor fllTf , 

We show that these apparently conflicting data might 
well be reconciled by a deeper understanding of how 
microstructure affects screening. It is well known, for 
example, that the screening length in weakly coupled 
Josephson array of grains is dominated by the magni- 
tude and temperature dependence of the Josephson cou- 
pling current between array elements fl48| . Thus, tunnel 
coupling across grain boundaries and/or planar defects 
(weak links), rather than the BCS response of the grains 
themselves, mainly determines the magnetic screening 
length, surface resistance, and critical current (see a re- 
view article Jl9]]). The extrinsic effect due to the weak 
links can lead to a linear T dependence in the effective 
A a b(T) at low temperatures and to a large residual sur- 
face resistance |30|. Similarly, Hebard et al. showed 
that the current-induced nucleation of vortex-antivortex 
pairs at defects can make an additional extrinsic contri- 
bution to the screening length, i.e., a pinning penetration 
depth X p ab (T). Within this scenario, the X v ab (T) in zero 
magnetic field is given by |0| 



K b (t) = Kb(o)/a-n 



(9) 



where t = T/T c ; X p ab {0) = [® / H c (0)}^/2NdF; $ is the 
flux quantum; Nd is the areal density of uniformly dis- 
tributed defects; H c (0) is the zero-temperature critical 
field. In the presence of the external dc field H, the ex- 
pression for X^ b (0,H) has to be modified |5(J. The total 

screening length is X a b(t) 



,TOi)] 2 + [A^(*)] 2 , where 
X^ b (t) is the intrinsic London penetration depth ^l]. As- 
suming an s-wave pairing symmetry, we readily show 
that the X ab (T) at low temperatures (below 0.2A(0)//cs) 
is given by 



Aah(T) = A a6 (0) + 
exp[-A(0)/A: B r] 



[A^(0)] 2 



y/nA(0)/2k B T 



Aab(0) 

, A 2 b (0) - [A^(0)] 2 
A ab (0)T c 2 



(10) 



It is clear that the T 2 dependence of X ab (T) at low 
temperatures in zero field can be completely caused by 
the extrinsic effect, that is, the nucleation of vortex- 
antivortex pairs at defects. If is negligible, X ab (0) = 
A^ b (0), and the second term in Eq. is absent. Then 
we recover the BCS expression |p2[, 



Kb(T) = X ab {0) + A ab (0) \ArA(0)/2fc B T exp[- A(0)/fc s T] . 

(11) 

In Fig. ^a, we plot temperature dependence of X a b(T) 
below 6 K for a PCCO single crystal (the data are from 
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Rcf. |47]]). The zero-temperature in-plane penetration 
depth A ah (0) was measured to be 2500 A @. This 
crystal shows T° nset at 22 K (defined by the onset of 
diamagnetism), and T"" d at 19 K (defined as the in- 
flection point on X a b(T)) A wide superconducting 
transition in this crystal manifests a rather low quality 
of the crystal. 

We fit the data by Eq. |l^ with two fitting parameters 
A(0) and A^(0), and with a fixed T c = 20.5 K (the aver- 
age of T° nset and T c mid ). The solid line is the fitted curve 
by Eq. [l0[ It is remarkable that the fit is very good. This 
can be seen more clearly in Fig. where the difference 
between the data and the fitted curve is plotted. There 
is a negligible systematic error (the deviation is less than 
the magnitude of the data scattering). From the fit, we 
find A{0)/k B = 29.6±0.1 K, and A^ b (0) = 1643 A. The 
deduced A^ b (0) is in excellent agreement with the value 
(1600±100 A) obtained from the optical data Q. The 
magnitude of 2A(0) /fcsT c = 2.9 is also the same as the 
one deduced from a high-quality film where the T 2 term is 
absent (46). The value of A(0) justifies the fit to the data 
below 6 K, namely, 0.2A(0)/fc B . Therefore, the X ab (T) 
data for the crystal are in quantitative agreement with 
an anisotropic s-wave pairing symmetry with no nodes. 
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FIG. 7. (a) Temperature dependence of X a b(T) below 6 K 
for a PCCO single crystal. The solid line is the fitted curve 
by Eq. [to]with 2A(0)/fc s T c = 2.9 and A^(0) = 1643 A. The 
value of A^(0) was found to be 1600±100 A from the optical 
data [^3|. (b) The difference between the data and the fitted 
curve. The data are from Ref. ^7j. 

From the values of A^ b (0) and A a fc(0), we calcu- 



late A^ h (0) = 1884 A. Using the relation A^ b (0) = 
[<f>o/H c {0)]y/2N d /n and H c (0) = 2 kOe @, we estimate 
N c i = 5.2 x 10 10 /cm 2 , corresponding to one defect over 
1333 Cu sites. This implies that a small density of de- 
fects can produce a quite large A^ h (0) which contributes 
a substantial T 2 term in A a f,(T). 
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FIG. 8. The T 2 dependence of 

the quantity l-A^(0)/A^ b (T) over 0.4-10.8 K for the same 
PCCO crystal as the one in Fig. ^. The crossover from the 
T 2 to a higher power-law dependence starts at about 5 K. 
There is no crossover from the T 2 to the T dependence at T* 
~ 9 K. 

In order to rule out the possibility that the data 
can be also consistent with a d-wave symmetry in the 
dirty limit, we plot the data as 1 — X 2 ab (0) / \ 2 ab (T) vs 
T 2 in Fig. ^. It is apparent that the quantity 1 — 
A 2 b (0)/A 2 6 (T) is proportional to T 2 below about 5 K. 
For a dirty <i-wave superconductor, a crossover from 
T 2 to T dependence should be seen at a temperature 
T* ~ A a6 (0)ln2/[A M (0)dA afc /dT 2 ], where A M (0) is the 
maximum gap at zero temperature f35f . Using A a b(0) 
= 2500 A M , d\ ab /dT 2 = 3.7 A/K 2 ||47), and A M (0) 
= 2.5T C |5j], one finds T* ~ 9 K. There is no such a 
crossover at any temperatures up to 11 K (see Fig. ||). 
Only a possible crossover from the T 2 to a higher power- 
law dependence is seen at about 5 K. Therefore, the 
data cannot agree with the ci-wave pairing symmetry. 
Furthermore, the absence of the linear T term in X a b(T) 
indicates that the extrinsic contribution to X a b(T) due to 
weak links (mJ is negligible in this crystal. 

In Fig. we show [A Qb (T) - A afc (0)]/A a6 (0) as a func- 
tion of temperature for a high-quality PCCO thin film 
(the data are from Ref. (46) ). The film has the lowest 
residual resistivity (< 50 fifl cm) and the highest T c (24 
K) reported for the PCCO system jh|. This indicates a 
high-quality of the film, which was grown using molecu- 
lar beam epitaxy. The optimal quality of the film may 
be due to the fact that a homogeneous oxygen reduction 
can be easily achieved in thin films. Since the data at 
low temperatures are quite flat, it appears that there is 
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neither T 2 nor T contribution. We thus fit the data be- 
low 6.5 K by Eq. [ll] with one fitting parameter A(0). 
The best bit gives A{0)/k B = 31.9±0.1 K, which justi- 
fies the fit to the data below 6.5 K (~ 0.2A(0)/fc B ). This 
leads to 2A(0)/fcsT c = 2.7, which is nearly the same 
as that deduced above for the less ideal crystal where 
there is a significant T 2 term in A a &(T) due to the ex- 
istence of defects. All these results consistently suggest 
that the pairing symmetry in electron-doped cuprates is 
the anisotropic s-wave with no line nodes. 
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FIG. 9. Temperature dependence of 

[X ab (T) - A o &(0)]/A o &(0) for a high-quality PCCO thin film 
with the lowest residual resistivity and the highest T c . The 
solid line is the fitted curve by Eq. |ll| with 2A(0)/fc B T c = 
2.7. The data are from Ref. Efl|. 

Polarized Raman scattering has also shown 

that the symmetry of the order parameter in 
Ndi.84Ceo.i6Cu04_j, (NCCO) is consistent with an 
anisotropic s-wave. More precisely, the tunneling spec- 
tra |3?j are consistent with a gap function: A(0) = 
A(s + cos40) with s > 1. If we use A M = 2.5T C Q and 
the minimum gap A m = 1.4T C (from the X a b(T) data), we 
find A(0) = 1.15(3.52+cos40) meVfor an electron-doped 
cuprate with T c — 24 K. Therefore, three bulk-sensitive 
experiments consistently indicate an anisotropic s-wave 
pairing symmetry in electron-doped cuprates. 

IV. PHASE-SENSITIVE EXPERIMENTS ALONG 
THE C-AXIS DIRECTION 

The most reliable phase-sensitive experiment is the 
atomically clean BSCCO Josephson junctions between 
identical single crystal cleaves stacked and twisted an an- 
gle 4>o about the c-axis Q. The quality of the junction 
is nearly the same as that of the intrinsic Josephson junc- 
tions made of the Bi-0 insulating layers. Theoretically, 
it has been shown that the critical current I c of the twist 
junction is |p8| 



I c cx ^2 Vi^i cos^o, (12) 

l 

where £ = 0, 1, 2, , and <q t « rjo for £>l. The 

above equation indicates that the s-wave component con- 
tributes to the critical current much more effectively. The 
experiment shows Jl^ j that the I c value is nearly indepen- 
dent of the twist angle <pQ, and the temperature depen- 
dence of I c is consistent with the Ambegaokar-Baratoff 
(AB) model for an s-wave superconductor. This indi- 
cates that the s-wave component in this material must 
be significant compared with the other high angular mo- 
mentum components. For slightly overdoped BSCCO, we 
have found that the gap function is A(0) = 17. 75 (cos 40+ 
0.46) meV for T c = 90 K, and A(0) = 18(cos40 + 0.33) 
meV for T c = 86 K. Then we have A s = 6-8 meV, which is 
not small compared with the g-wave component A g = 18 
meV. Since 774 << 770 |58[| , the dominant contribution to 
the I c should be the s-wave component, as observed [ p^[ . 
From the magnitude of the s-wave component, we can 
calculate I c Rn = (7r/2e)A s = 9-12 mV. The measured 
I c Rn value is about 8 mV p2| . This is in quantitative 
agreement with the predicted value considering the fact 
that the strong coupling effect can reduce the I c Rn value 
by more than 20%. 

Another reliable phase-sensitive experiment is the c- 
axis Pb/YBa2Cu307_j / Josephson junction experiment 
p3[ . This junction can be described as SINS' (where S 
= YBCO, S'= Pb, I and N represent the insulating and 
normal- metal layers, respectively). Due to a very short 
coherent length £ c along the c-axis direction, the bulk gap 
will be strongly depressed at the SI interface; the depres- 
sion factor is c/£ c (where c is the lattice constant along 
the c-axis) From £ c = £, a b/l (where 7 is the mass 

anisotropy parameter and equal to about 8 for optimally- 
doped YBCO @), we get £c = 1-7 A by taking £ afc = 
14 A. Therefore, the gap size at SI interface will be sup- 
pressed by a factor of about 7. Since the bulk s-wave 
component A s in slightly overdoped YBCO is 3-5 meV 
(see above) , this component at the SI interface should be 
reduced to 0.4-0.7 meV. Then the I c Rn value is calcu- 
lated to be 0.93-1.27 mV, in quantitative agreement with 
the measured one (~0.9 mV) p^ |. 

Now we discuss another c-axis Josephson tunneling ex- 
periments in which a conventional superconductor (Pb) is 
deposited across a single twin boundary of a YBCO crys- 
tal |Q. Because Pb is an s-wave superconductor, the Pb 
counterelectrode couples only to the s-wave component 
of the YBCO order parameter. If YBCO were predomi- 
nantly G?-wave, any small s-wave component added to the 
dominant d-wave component would change sign across 
the twin boundary. In this case, magnetic fields parallel 
to the boundary would produce a local minimum in I c 
at B = 0, in agreement with the observation f6l)| . The 
experimental results thus appear to provide evidence for 
mixed d- and s-wave pairing symmetry in the bulk with 
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a reversal in the sign of the s-wave component across the 
boundary. However, if the bulk OP symmetry in a single 
domain were d+s or d—s, one would expect a nearly zero 
I c in heavily twinned crystals. The fact that the observed 
I c Rn in heavily twinned crystals |]l3f is nearly the same 
as the one in the single-domain crystal rules out the 
bulk d + s or d — s wave OP symmetry in YBCO. There- 
fore the only possibility is that a half/fractional flux is 
trapped in the twin boundary. Also, this can naturally 
explain why I c does not go to zero even for a symmetric 
junction with the same junction area in both sides of the 
twin boundary J6l[ |. 

V. PHASE-SENSITIVE EXPERIMENTS ALONG 
THE AB-PLANES 

The phase-sensitive tricrystal experiments on both 
hole and electron-doped cuprates |2^,^,|l^| show that the 
OP symmetry is the d-wave, in contradiction with the 
above conclusion drawn from many bulk-sensitive exper- 
iments. In order to resolve the above discrepancy, one 
should notice that the tricrystal experiments are rather 
surface sensitive, so these experiments are probing the 
OP symmetry at the surface/interface, rather than in 
the bulk. Based on the Ginzburg-Landau free energy, 
Bahcall [jl8| has shown that the OP symmetry near sur- 
faces/interfaces can be different from that in the bulk 
if the bulk OP is strongly suppressed at the surfaces. 
Experimentally, the observed I c Rn values in all the 
tricrystal experiments are about two order of magnitude 
smaller than the intrinsic Ambegaokar-Baratoff limit. 
For example, in the optimally-doped YBa2Cu307_j,, the 
magnitude of the maximum gap Am(0) is about 30 meV 
Then the intrinsic I c Rn value should be equal 
to the Ambegaokar-Baratoff limit 7rAM(0)/2e = 47 mV, 
which has been confirmed by a nearly ideal SIS break 
junction experiment fl7}| . However, the observed I c Rn 
values in the tricrystal experiments on YBa2Cu307_j, 
and Tl 2 Ba 2 Cu0 6+ .y §,[l0) are about 1.8 mV and 0.5 mV, 
respectively. These values are about two order of magni- 
tude smaller than the intrinsic bulk values. Similarly, the 
observed I c Rn value in the NCCO and PCCO tricrystal 
experiments is about 0.1 mV, as inferred from the mea- 
sured critical current density J c = 6 A/cm 2 |^] and the 
empirical relation between I c Rn and J c fl62| . This I c Rn 
value is also about two order of magnitude smaller than 
the intrinsic bulk value, which is estimated to be ~ 8 mV 
with A M (0) = 2.5T C @. Therefore, the OP at the inter- 
faces of the grain boundary junctions must be strongly 
depressed in order to explain such small I c Rn values. 
This strong depression in the order parameter ensures 
the condition under which the OP symmetry near sur- 
faces/interfaces can be different from that in the bulk 
p8[ . Hence, it is very likely that the tricrystal experi- 
ments are detecting the OP symmetry at the degraded 



interfaces, which may be different from the intrinsic one 
in the bulk. 

Now the question arises: why the bulk OP can be so 
strongly depressed at the surfaces of the grain bound- 
ary junctions? It is known that the coherent length 
in cuprates is very short due to a large superconduct- 
ing gap and small Fermi velocity. The short coherent 
length in cuprates can lead to a large depression of the 
OP near the interfaces even within the conventional the- 
ory of the proximity effect |^,^3|. Alternatively, several 
groups p4|-|66| showed that there are possibly nonsuper- 
conducting regions near the boundary of the junction due 
to hole depletion and/or strain, so that the critical cur- 
rent density can be reduced by several order of magnitude 
compared with the intrinsic bulk value. 

There is another way to explain the tricrystal experi- 
ments. As discussed above, the bounaries of the grain- 
boundary junctions are intrinsically underdoped super- 
conductors or nonsuperconducting antifcrromagnets due 
to hole depletion and/or strain |64|-|67|. For underdoped 
cuprates, the superconductivity mainly arises from the 
Bose-Einstein condensation of preformed pairs pq] . In 
this case, the symmetry of the superconducting conden- 
sate is different from the pairing symmetry; the former is 
the e?-wave while the latter might be s-wave ]6q| . Since 
Josephson tunneling probes the symmetry of the super- 
conducting condensate, the d-wave symmetry of the con- 
densate is consistent with the tricrystal experiments. 

VI. CONCLUSION 

In conclusion, the existing bulk and nearly bulk- 
sensitive experiments consistently show that the domi- 
nant bulk pairing symmetry in hole-doped cuprates is of 
extended s-wave with eight line nodes, and of anisotropic 
s-wave in electron-doped cuprates. The deduced ex- 
tended s-wave pairing symmetry for hole-doped cuprates 
is also in quantitative agreement with the phase-sensitive 
experiments along the c-axis for both Bi2Sr2CaCu208+ 2/ 
and YBa2Cu307_ y . The proposed pairing symmetries do 
not contradict some surface- and phase-sensitive experi- 
ments which show a predominant d-wave pairing symme- 
try at the degraded surfaces. 
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